Orbital evolution of a test particle around a black hole: 
Indirect determination of the self force in the post Newtonian approximation 
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Comparing the corrections to Kepler's law with orbital evolution under a self force, we extract the 
finite, already regularized part of the latter in a specific gauge. We apply this method to a quasi- 
circular orbit around a Schwarzschild black hole of an extreme mass ratio binary, and determine 
the first- and second-order conservative gravitational self force in a post Newtonian expansion. We 
use these results in the construction of the gravitational waveform, and revisit the question of the 
relative contribution of the self force and spin-orbit coupling. 
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I. INTRODUCTION 

The problem of finding the gravitational self force for 
an extreme mass ratio binary has drawn much atten- 
tion over the last decade, especially because of its impor- 
tance for the calculation of gravitational waveforms. The 
usual approach for the calculation of the self force is to 
calculate it directly, via a number of currently available 
methods In any direct calculation, one of the major 
difficulties is the regularization problem of the self force, 
which is currently not under full control although over 
the last few years major new insights have been achieved 
towards its solution. Finding the self force for generic or- 
bits requires such direct calculations. Once the self force 
is found, one needs to find the resulting orbital evolution, 
and to find the waveforms. Direct orbital evolution (for 
quasi-circuar orbits in Schwarzschild, for nonspinning |2j 
and spinning particles [2j) has been done with a scalar 
field self force, but hasn't been done yet for the gravita- 
tional self force. 

The most immediate problem for evolving the orbit 
under the gravitational self force is that we do not have 
as yet the latter calculated via the direct approach. An 
alternative approach for finding the self force is an indi- 
rect one. Specifically, one can use quantities whose effect 
on the orbital evolution does not require regularization — 
or whose regularization is already understood — and ex- 
tract from them the finite, already regularized, part of the 
self force. A simple application of this approach is well 
known: One can start with the flux of energy and angular 
momentum radiated to infinity (and down the horizon of 
the central black hole), use a balance argument to find the 
rate of change of the particle's constants of motion, and 
use the latter to find the dissipative part of the self force 
0. This indirect determination of dissipative self force 
effects requires no regularization. As we show in this 
paper, extracting the conservative piece of the self force 
indirectly does require regularization, but none beyond 
the regularization that is needed for the derivation of 
the quantities that are used for the extraction of the self 
force. However, as is well known, the conservative part of 
the self force, which is of much practical importance for 



the accurate determination of gravitational waveforms, 
cannot be found using the flux of fields. We therefore 
extract it using other quantities, specifically the angular 
velocity. 

As the self force we derive is extracted from the an- 
gular momentum, there is no new information in it be- 
yond that that already is encoded in the angular mo- 
mentum. Therefore, the orbital evolution (and the wave- 
forms) obtained here are not different from those ob- 
tained by other methods that are based on the angular 
momentum. Specifically, in this paper we make use of a 
post-Newtonian expression for the angular momentum. 
Consequently, our waveforms are no different from post- 
Newtonian waveforms. The new content of this paper 
is in the viewpoint: the orbital evolution is obtained by 
a direct integration of the gravitational self force. To 
do this, we demonstrate how the (conservative piece of 
the) self force can be found from the angular momentum 
expression. 

In this paper, we consider the determination of the 
conservative piece of the self force using an indirect 
method — which in that sense is analogous to the determi- 
nation of dissipative self forces using fluxes — and apply it 
to a simple situation, specifically that of a quasi-circular 
planar orbit around a Schwarzschild black hole, i.e., or- 
bits that would be circular in the absence of radiation re- 
action. Specifically, we make use of (the post Newtonian 
correction of) Kepler's law to find not just the first order 
conservative self force, but also the second order conser- 
vative self force. After finding the self force in a specific 
gauge, we can use it to find orbit integrated gauge invari- 
ant quantities, such as the waveform. These results can 
be used for testing and checking direct results for the con- 
servative piece of self forces, and the corresponding reg- 
ularization schemes (including second order ones). We 
emphasize that the waveforms themselves are not new, 
as they are equivalent to the waveforms obtained by a 
post-Newtonian approach that keeps conservative effects 
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This paper follows previous papers on the orbital evo- 
lution of extreme mass ratio binaries under conservative 
self forces. In Ref. 0] — henceforth paper I — we showed, 
using quasi-circular orbits in Schwarzschild, that conser- 
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vative self forces are required for finding the waveform to 
O(e ), where e is the mass ratio of the binary. It was also 
shown that second order (dissipative) self forces, i.e., self 
forces to 0(e 3 ), are required for a full, self consistent de- 
termination of the waveform to this order. In Ref. 
hereafter paper II — the question of the relative impor- 
tance of conservative forces and spin-orbit coupling was 
addressed, and in Ref. |6( the equations of motion were 
solved perturbatively also for spin-spin coupling, again 
for quasi-circular orbits. In this paper we show how one 
can find the conservative self force, and the resulting or- 
bital evolution and waveform, from the post Newtonian 
equations of motion, not only for the first-order self force 
(i.e., self force to (3(e 2 )), but also the second-order self 
force (i.e., self force to 0(e 3 )). As pointed out in paper I, 
the latter are of importance for a full description of the 
waveform, and have gained recently much attention 0. 
The method of this paper, however, is not limited to the 
post Newtonian approach. 

The organization of this paper is as follows: In Section 
II we recapitulate the perturbative solution of paper I for 
the equations of motion, and extend it to the required or- 
der. In Section III we use this solution to find the self 
force indirectly. In Section IV we use this self force to 
find the orbital evolution and the waveforms, and finally, 
in Section V we revisit the question of the relative im- 
portance of self forces and spin-orbit coupling. 



II. PERTURBATIVE SOLUTION FOR THE 
EQUATIONS OF MOTION 

Let a non-spinning test body of mass fj, be acceler- 
ated because of its self force in a quasi-circular or- 
bit around a Schwarzschild black hole of mass M (with 
e := fi/M <C 1), such that its four- velocity is u a and its 
proper time is r. The equations of motion (EOM) are 
given by 



Du a 
~dV 



SF a/3 



J (3 9 



(1) 



which we solve perturbatively. 

The metric in the usual Schwarzschild coordinates is 
given by 



-F(r) dt 2 + F- 1 (r) dr 2 + r 2 dfl 2 



(2) 



where F(r) = 1 - 2M/r. Here, dti 2 = d-d 2 + sin 2 fid(j) 2 . 

As the orbit is planar, the 9 component of the EOM 
is trivial. We use the normalization condition for u a , 
namely u a u a = —1, to eliminate u l from the EOM 
Du l I dr = /i _1 /f F g lfc , i,k = t,r,4>, where D denotes 
covariant differentiation compatible with the metric J5J). 
We next use the t component of the EOM to eliminate 
m*. We can simplify the EOM to first-order (nonlinear) 
ODEs by taking r = V(r), x = Vx'(r), where x denotes 
any quantity. We find the EOM to be 
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where {u 1 ) 2 = 1/[1 - Ml/r - r 2 a - V 2 /(l - 2M/r)]. We 
denote by the orbital frequency, and by an overdot and 
a prime (partial) differentiation with respect to coordi- 
nate time t and r, respectively. Here, a measures the 
deviation from Kepler's law, i.e., 



VL 2 = M/r 3 + cr(r) 



(5) 



This last relation is gauge independent, although each 
of the terms on the right hand side are separately gauge 
dependent. Specifically, we fix the gauge by choosing the 
first term on the rhs to include the radial Schwarzschild 
coordinate r: under an infinitesimal coordinate trans- 
x' a = x a + £ a where the gauge vec- 



h! 
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formation x a — > x' a = a 
tor £ Q = 0(ji), the metric perturbations h a p 
h a p — 2V( Q £ / 3). In the last expression, the covariant 
derivative operator V Q is compatible with the back- 
ground metric @- A gauge choice on the metric pertur- 
bations h aj 3 (e.g., the Regge- Wheeler or Lorenz gauges) 
is therefore equivalent to a condition on the gauge vec- 
tor The latter may have a radial component, that 
changes the radial coordinate description of the orbit. 
As noted by Detweiler ,8], while the angular velocity fl 
is gauge invariant, the radius of the orbit r is not. This 
means that under a gauge transformation the different 
terms on the rhs of Eq. change, but in such a man- 
ner that their combination, or the lhs of is invari- 
ant. One may, therefore, fix the gauge (possibly up to 
a residual gauge freedom) by determining the ratio of 
the terms on the rhs of JSJ, or, alternatively, by fixing 
the geometrical meaning of the symbol r appearing on 
the rhs of (JSJ. In particular, the latter may be fixed 
to equal the Schwarzschild radial coordinate of the un- 
perturbed Schwarzschild background (0). By clarifying 
the geometrical meaning of r we therefore effectively fix 
the gauge. Notice that this fixing is not equivalent to 
merely choosing which coordinates are used to describe 
the background geometry. 

We next expand in powers of e as a(r) = + a^ 2 \ 
V = + V^ 2 \ and a t = af ) + af\ x^ denoting 
the term in x which is at 0(e J ), and being the (self) 
acceleration. We then expand the self force as /• = 

ft + , where f^ = ^af \ Solving perturbatively, 
we find the first order terms to be 
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and the second-order terms to be 
r-3M 



(7) 



r(2) 



(2) 



(r - 2M) 2 



r mm 3MFW 2 v(Dv(iy 

/i ,/r r(r - 2M) 2 r - 2M 



(8) 



F (2) = r(r - 3M) 
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In papers I and II we were primarily interested in the 
orbital evolution, and therefore attempted to justify ne- 
glecting the (dissipative) terms in Eq. © that are linear 
in the second-order self force. In this paper we continue 
to assume this to be the case. Unlike the prequels, in this 
paper we introduce Eq. (JSJ, and leave there the (conser- 
vative) term linear in the second-order self force. In fact, 
we make use of Eq. |(HJ| to isolate this term, and find in- 
directly the (post Newtonian) expression for the second- 
order (r component of the) self force. The first-order 
counterpart is found similarly from Eq. ©. 



III. FINDING THE (POST NEWTONIAN) SELF 
FORCE INDIRECTLY 



Consider the correction to Kepler's third law, in a post 
Newtonian expansion: 



l2 _ M + p 



n z = 



i? 3 



I f ! -■'-! - /'h - j 6 + ^v + v 2 ) 7 2 
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+ ( -10 + pv+ —v z + v A ) Y 



(10) 



where R is the harmonic radial coordinate, v := 
Mp/(M + p) 2 : 7 := (M + p)/R, and p is a certain 
(known) expression |9J. Note, that while Vi 2 is a gauge 
invariant quantity, it is here expressed in terms of gauge 
dependent quantities, specifically the harmonic coordi- 
nate R. 

This correction of Kepler's law is conservative. Equa- 
tion (|10|) was derived in [jj using Hadamard Partie finie 
regularization, and is correct for arbitrary mass ratios, 
although we make use of it here for v <C 1. We further 



comment that as the self force we extract below is ob- 
tained from Eq. I|10|) . there is no new physics in the self 
force or the resulting orbital evolution beyond that that 
is already in Eq. (|10[1. 

We first recognize that the t'-independent terms inside 
the square brackets in Eq. (|10|) are the leading terms in 
the expansion of (1 + 7) -3 for 7 C 1. Using this, we 
notice that i?~ 3 (l +7)~ 3 = r~ 3 [l - 3p/r + 0(p 2 )}. With 
this substitution, we express the angular velocity as 



M 



M/i 61 Af 2 u 
65 Mp 2 



(11) 



where following the Keplerian term on the rhs we present 
in the first line of 1|11|) the Newtonian self force correction 
[Tfl |. followed by the 1PN and 2PN first-order self- force 
corrections, and where in the second line we present the 
Newtonian and 1PN second-order self-force corrections. 

The correction to Kepler's law can be construed as 
caused by a conservative self force, that we decompose in 
the post Newtonian framework to the K th PN order as 



U) 



K 

E 

fe=0 



A3) 

J% kPN ■ 



(12) 



In this decomposition we formally treat the Newtonian 
self force as a zeroth order post Newtonian term. 



Substituting to 2PN crW 



p/r 



(61/4:)M 2 p/r 5 in Eq. ©, we find that 
■ 2 p 2 M 73 p 2 M 2 



2Mp/r 4 



+ 0{p 2 M i /r b ). (13) 



This is our starting point in the investigation of the or- 
bital evolution under the gravitational self force in the 
next section. 

Next, we find the radial component of the second- 
order self force. Substituting to 1PN = -3p 2 /r 4 + 
(65/4)Mp 2 /r 5 in Eq. ©, we find that 



37 p 3 M 



+ 0(p 3 M 2 /r 5 ). 



(14) 



To this PN order, the dissipative terms in the rhs of 
Eq. ((SJ do not contribute, as they enter at 2.5PN order. 
We can therefore ignore the terms involving dissipation, 
i.e., all the terms including U^ or /j in the extraction 
of^ 2) . 

We note that the result for f r , although derived from a 
gauge invariant quantity, certainly depends on the choice 
of gauge. For example, it was derived in Ref. |ll| in the 
Regge- Wheeler gauge to 1PN. Notably, even the ques- 
tion of what is the sign of f r (and, consequently, also 
the question of whether it is attractive or repulsive) is 
gauge dependent. The only meaningful quantities to find, 
therefore, are gauge independent ones, such as f2 2 . In- 
deed, when the expression of Ref. [ll| is used to find the 
angular velocity, Eq. i|ll[l is recovered in the appropriate 
PN order. 
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IV. ORBITAL EVOLUTION AND WAVEFORMS 

Armed with Eq. (|13f) we can now integrate the orbit 
and find the gravitational waveforms. We use here the 
method of paper II with only slight modifications, and 
the reader is referred there for detail. The dissipative 
part of the self force is modeled at large distanced by the 
3.5PN flux of energy to infinity [^j, and at small dis- 
tances by the numerically computed luminosity in grav- 
itational waves |l3j| , which we fit to a polynomial to find 
a smooth function, from which the derivatives of can 
be found. As the relative error in the determination of 
the luminosity in either method becomes comparable at 
r/M ~ 18, we make a transition from one to the other 
at that distance. Specifically, we use a Fermi-Dirac type 
mixing, of the form 



dE 11 



dt 



p(r) 



dE 



PL 



dt 



[l-p(r)] 



dE 3 



.5PN 



dt 



(15) 



where p(r) = l/{exp[(r — ro)/A] + 1}. In practice, we 
took r = 18M and A = 10~ 3 M. We find this type of 
mixing to give better results than the sliding average we 
used in paper II. 
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FIG. 1: The waveform for an orbit starting at r — 10M at 
t = 0, and decaying to the ISCO at 6M for fi = 10~ 4 M. 
Dotted line: no self force. Dash-dotted line: Newtonian self 
force. Dashed line: 1PN self force. Solid line: 2PN self force. 

As in papers I and II, the orbit is found through 



t(r) 



dr 



and ip(r) 



O(f) dr 



(16) 

and the waveform is obtained through the usual "re- 
stricted waveform" approximation (see paper I for more 
detail). Figure Q] shows the waveforms for an orbit that 
starts at r — 10AT, and decays to the inner-most stable 
circular orbit (ISCO) at r — 6M for self-force correction 
at the Newtonian, 1PN and 2PN orders. Not surpris- 



ingly, close to the ISCO the post Newtonian approxima- 
tion breaks down. 
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FIG. 2: The waveform for an orbit starting at r — 10M at 
t = 0, and decaying to the ISCO at 6M for fi = 10" 4 M 
at 1.5PN order for the conservative effects. Solid line: non- 
spinning body. Dotted line: Maximal spin aligned with the 
orbital angular momentum. Dashed line: Maximal spin anti- 
aligned with the orbital angular momentum. 



V. ON THE RELATIVE CONTRIBUTIONS OF 
SELF FORCES AND SPIN FORCES 

In paper II we considered the relative importance for 
orbital evolution of the conservative self force and effec- 
tive forces that come about because of spin-orbit cou- 
pling. The argument of paper II was that for a spinning 
body the spin-orbit effective force is quadratic in the 
body's mass, and therefore should be included in the con- 
struction of waveforms for any self-consistent description 
of the latter at O(e ). Specifically, paper II considered 
a conservative self force that is modeled after the scalar 
field self force that was found in [14( . Paper II then sug- 
gested that under this assumption, the self force effects 
on the waveforms are likely to be overwhelmed by the 
spin-orbit coupling effects. 

More recently, Pound et al have correctly commented, 
that the model adopted in paper II was irrelevant for 
realistic binaries, because the gravitational self force en- 
ters at a lower PN order than the scalar field counterpart 
(Tlj . Indeed, the conservative gravitational self force is a 
1PN effect (not counting here the Newtonian self force), 
whereas the conservative scalar field self force is a 3PN 
effect 14] . On the other hand, the spin-orbit force is 
a 1.5PN effect. It might be inferred from the argument 
by Pound et al that for realistic binaries the situation is 
reversed: it is not the self force that can be neglected 
compared with the spin-orbit effect, but the other way 
around. 
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In this Section we revisit the question of the rela- 
tive importance of the conservative self force and spin- 
coupling effects. We follow the method of paper II, and 
include the spin-orbit force in the EOM. Figure [21 dis- 
plays the waveform at 1.5PN order for the conservative 
effects for a non-spining body, and bodies with maximal 
spins aligned or anti-aligned with the orbital angular mo- 
mentum. Bodies with intermediate values for their spin 
will be described by waveforms in between these two. 
Comparison of Fig. and [21 suggests that the orbit in- 
tegrated spin-orbit effect is comparable to the self force 



effect, even though the former enters at a higher PN order 
than the latter. We therefore infer that any self consis- 
tent determination of the waveform to O(e ) will have to 
include both the self force and the spin-orbit effect. 
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